How Metallic are Small Sodium Clusters? 
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Cryogenic cluster beam experiments have provided crucial insights into the evolution of the metal- 
lic state from the atom to the bulk. Surprisingly, one of the most fundamental metallic properties, 
the ability of a metal to efficiently screen electric fields, is still poorly understood in small clusters. 
Theory has predicted that many small Na clusters are unable to screen charge inhomogeneities and 
thus have permanent dipole moments. High precision electric deflection experiments on cryogeni- 
cally cooled Najv (N < 200) clusters show that the electric dipole moments are at least an order 
of magnitude smaller than predicted, and are consistent with zero, as expected for a metal. The 
polarizabilities of Na clusters also show metal sphere behavior, with fine size oscillations caused by 
the shell structure. 



By definition, a classical metal is a material which can- 
not support an internal electric field. An electric field 
E ext (r) applied to a metal object of arbitrary shape will 
cause the charge density to rearrange so that E; nt = 0. 
A caveat of this property is that a metallic object cannot 
have a permanent electric dipole moment (or any other 
moment), since this implies that there is a non- vanishing 
internal electric field [T] . This property of metals applies 
on the macroscopic level, but it is not a priori obvious 
that it applies to extremely small objects such as metal 
clusters. The effectiveness of the screening can be experi- 
mentally tested by measuring the electric dipole moments 
and polarizabilities. 

Early experimental and theoretical work on metal 
clusters focused on the static dipole polarizability and 
demonstrated that alkali metal clusters could be approx- 
imately treated as small metal spheres [2] . This led to the 
well-known jellium model which allowed a self-consistent 
description of the electronic shell structure of small clus- 
ters . The spherical jellium model predicts that the po- 
larizability of an alkali cluster is a(N) = (R + S(N)) 3 , 
where R — r^A 1 / 3 is the classical cluster radius, r s is 
the Wigner-Seitz radius, N is the cluster size (in atoms). 
8{N) is a quantum correction to the radius, often re- 
ferred to as the spillout factor since it indicates that the 
electronic screening actually extends beyond the classi- 
cal cluster radius. To first order, 8(N) is constant and 
comparable to the Lang-Kohn value for jellium surfaces 
[3]. In more sophisticated calculations, S(N) varies with 
cluster size and shows non-trivial shell structure effects 
0H6]. 

The spherical jellium model is clearly flawed: a small 
metal cluster is not even approximately spherical [2] [7] , 
and the ionic structure has been shown to have signifi- 
cant effects on the thermodynamic properties, [51 [S] and 
photoelectron spectra [10] . Nevertheless, many physi- 
cal properties, including the polarizabilities |11H13| are 
surprisingly well-described. The existing experimental 
data on Na cluster polarizabilities only sparsely covers 
the range of cluster sizes, and the experiments were done 
at temperatures where the clusters are liquid [SJ [5] • How- 



ever, as we show here, essential features of the jellium 
model are still observed even in high precision measure- 
ments, at cryogenic temperatures where the clusters are 
expected to be rigid with few or no excited vibrations (20 
K). 

Electric dipole moments are also expected to be highly 
sensitive to the electronic screening, and dipole moments 
have been observed in many metal cluster systems (e.g. 
Nb, V, Ta Q31 H5], and Sn and Pb QH [TF] ) . An asym- 
metric cluster without inversion symmetry is expected 
to have an electric dipole moment, and its magnitude de- 
pends on how effectively the charge inhomogeneity of the 
ion cores is screened by the valence electrons. In the case 
of Pb clusters, the link between reduced screening and 
dipole moments is supported by a recent experiment [18] 
which found reduced core-hole screening in the same size 
range where dipole moments were observed [17j [19) . 
For small clusters this reduced screening has been ex- 
plained as a consequence of partial localization of the 
electrons due to change in the bonding character and 
low coordination[18]. An all-electron quantum chemi- 
cal calculation has predicted that similar effects will lead 
to dipole moments in Na clusters [20] . Our experiment 
shows that the electric dipole moments are much smaller, 
and that metallic screening is not well described by the- 
ory even for a cluster as small as Na3 . The failure of the- 
ory to correctly describe static screening in metal clusters 
is a serious outstanding problem. 

Electric dipole moments and polarizabilites are ide- 
ally measured using cryogenic molecular beam deflection 
methods. A beam of neutral metal clusters is produced, 
deflected and detected using methods that have been 
previously described [2] (see Refs. [141 ED f° r experi- 
mental details and parameters). Briefly, cryogenically 
cooled sodium clusters are produced in a laser vaporiza- 
tion (Nd:YAG 532 nm; 5 mJ/pulse) cluster source op- 
erating at 20 K. The beam velocity is measured with a 
mechanical chopper. The cluster beam is collimated ( 0.1 
mm slits) and passes through the pole faces of an inho- 
mogeneous electric field {E = 85 kV/cm, dE/dz=2l8 
kV/cm 2 ). The clusters then deflect due to the force 
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caused by the electric field gradient on the electric dipole 
that has an intrinsic component and an induced com- 
ponent. The induced dipole moment causes a uniform 
deflection of the cluster beam, while the intrinsic dipole 
moment (primarily) causes a broadening of the beam (see 
below for details). The cluster beam then enters a posi- 
tion sensitive time of flight mass spectrometer, that si- 
multaneously measures the mass and deflection for all 
clusters in the beam. This method has been previously 
used to measure the electric dipole moments of large po- 
lar molecules and clusters [2~T1. 
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FIG. 1. Electric deflection profiles for Najv N=2,6,18, and 
50. The bold line shows the beam profile with the field off, 
the light line is with the field on. The green dashed curve 
shows a simulation of the deflection profile using the adia- 
batic rotor model [22] , with the dipole moments calculated in 
Ref. [2D]. For Na6 a significant depletion of the beam inten- 
sity is observed which can be explained by an isomer with a 
large dipole moment. 

For every species in the beam, we measure a distribu- 
tion of polarizations. We assume that the induced polar- 
ization P = P a + P p is due to two effects: the electronic 
polarizability P a = aE, and the dipole moment p pro- 
jected onto the field, time-averaged over the rotational 
motion P p = (p z )t- For a metal spheroid, a is an aver- 
age of the principal polarizabilities [T2"] . Because P p de- 
pends on the initial conditions (orientation, energy, and 
angular momentum) when a cluster adiabatically enters 
the deflector, the ensemble of clusters shows a distribu- 
tion of polarizations p(P p )- Clusters will in general be 
deflected toward both the high and low field directions, 
depending on their initial orientation. The observed de- 
flection profile is thus a convolution of the beam profile 
with p(P), and the signature of the dipole moment is a 
broadening of the molecular beam, which we measure by 
Act = \fo% a — cr^ff. (where cr on / ff are the width of the 



peaks with the field on or off, respectively) 

To derive a quantitative relation between the dipole 
moment p and the beam broadening Act, we use the 
adiabatic rotor model developed by Bertsch and others 
[2"T1 |2"2"] . This model uses classical rigid-body mechanics 
to calculate P p = (p z )t- For Nai at 20 K, the rota- 
tional constant B = h 2 /2I w IfieV, so 2B/kT w 0.001, 
thus the rotational levels are effectively continuous, and 
classical mechanics applies. For a spherical rotor in the 
pE/kT <C 1 limit, the model predicts a polarization dis- 
tribution with the analytic form: p(P) = (l/2p) log \p/P\ 
[22] ■ The variance of this distribution is p 2 /9 so the de- 
flection profile of a cluster with p will show Ac = p/3. 
For our experiment, p = 0.1 D, gives pE/kT w 0.01, 
so the asymptotic regime pE/kT <C 1 applies. The 
structure of the cluster also effects the deflection pro- 
file. For symmetric tops (Ri ^ R2 = R3), the quantita- 
tive relation between p and Act is slightly different when 
pE/kT < 1. Simulations for R1/R3 = 1.4 show that 
the relation p = 3 Act holds to within 7%. Na clusters 
are known to show triaxial distortions, and there has 
been experimental and theoretical work |23) suggesting 
that a polar asymmetric rotors will tumble chaotically in 
the field if perturbed. This explanation was invoked to 
explain deflection experiments on biomolecules |23j with 
dipole moments of 6 D, that showed reduced broadening. 
In our laboratory, we have performed deflection experi- 
ments on weakly polar, highly asymmetric metal clusters 
(e.g. the planar Aug cluster (0.28 D) [24]) and observed 
no evidence of chaotic tumbling. In this case, the beam 
is still symmetrically broadened just as in the symmetric 
top case, and the p = 3Act estimate agrees with the value 
from multiple quantum chemical calculations. The model 
also assumes that any dipole moment is fixed in the clus- 
ters structure, and that the cluster is a rigid object. At 
20 K, the clusters are well below both the melting tem- 
perature and the range of temperatures where softening 
effects like premelting are known to occur [SJ E] ■ 

Per atom p/N and total p dipole moments estimated 
from the beam broadening using p = 3Act are shown in 
Fig. [2] Note that p/N scatters around zero for all clusters 
N > 20. For N < 20 there is a small amount of residual 
beam broadening which cannot be explained away as an 
artifact. [35] • For all of cluster sizes p/N is less than 
0.002 D per atom. 

The measured dipole moments appear to be greater 
than for N < 20 and Na3 appears to have the largest 
p/N. Yet, its total moment is only about 0.01 D. This 
measured value agrees with the measurements of Ernst 
[26] . and is significantly lower than Ref. [20 which pre- 
dicts a value of 0.3 D. It should be noted that there is 
agreement in the overall trend of the measured and calcu- 
lated dipole moments, which indicates that the calculated 
shapes could be accurate but that screening is severely 
underestimated, even for very small clusters. Despite the 
success of simple shell models, a Na cluster is a many- 
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the shell structure, which were not previously observed. 
(Fig. |3| Note that the minima in a/N correspond to 
spherical shell closings (e.g. lp 6 (N = 8), Id 10 (N = 18), 
l/ 14 (N = 34) , lg 18 (N = 58) , lh 22 {N = 92), 
and within measurement error, lh 22 (N — 186)). How- 
ever shell closings do not always correspond to minima 
in a/N. For example, maxima are observed for 2d 10 
(N = 68), 2/ 14 (N = 106), and perhaps 2g 18 (N = 156). 
Hence, the systematic trend is, shell closings with princi- 
pal quantum number 1 are minimum in a/N, and those 
with principal quantum number 2, tend to be maxima. 



FIG. 2. Experimental dipole moments for Na clusters at 20 K, 
estimated from the spherical rotor model p = 3-\/o"o n — cr^ s . 
(inset) Comparison between the theoretical dipole moments 
calculated in Ref. [2D] and the experimental values above for 
small Na clusters. Multiple theoretical values correspond to 
different isomers, dipole moments of zero (due to symmetry) 
were predicted for N=2,3,4,6,7,8,10, and 19. Thus the discrep- 
ancy is greatest for N = 12-19. The 1.2 D theoretical value 
for Na3 is a high energy isomer that should not be present in 
our 20 K experiment. 



body problem and practical calculations require approx- 
imations. Ref. [5D] uses a hybrid functional (B3LYP) 
to treat the exchange and correlation for the calcula- 
tion of a and p. This method is now known to show 
serious problems with bulk metallic systems |27j . There 
are many theoretical methods which deal with the many 
body problem at lower levels of approximation. Under- 
standing the origin of this error will require a comparison 
of these methods where the effect of the cluster struc- 
ture is carefully controlled for. Note that the electro- 
static energy of a Na2o cluster with a dipole moment 
of 0.1 D (far larger than what has been measured) is 
E = w 17/xeV. This suggests that correctly calcu- 
lating the charge density requires high energy accuracy. 

Further note that Ref [20 predicts two stable isomers 
for Nag. One is a planar triangle with p = 0, while the 
other is a pentagonal pyramid with p ss 0.5 D. Indeed, the 
observed intensity loss of Na6 (Fig. [I]) with applied field 
is consistent with two stable isomers, one which having a 
much larger p than the other. Experiments are planned 
to further investigate Na6- Note that for all other clusters 
there is no significant change in the total beam intensity 
when the electric field is turned on. 

We next turn to the high precision polarizability mea- 
surements for Na w 1 < N < 200. (Fig. [3}. First note 
that the measurements of a(N)/N generally agree with 
previous reports. The overall decreasing trend with in- 
creasing cluster size agrees with the simple approxima- 
tion for the polarizability of a conducting sphere with a 
spillout-cnhanced radius a(N) = (R + S) 3 . 

Besides the overall decreasing trend, the present mea- 
surement also clearly reveals variations in a/N with 
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FIG. 3. (upper) a/N of Najv (N = 1 - 30) at a beam tem- 
perature of 20 K. Compared with previous higher tempera- 
ture experiments [11H13] a/N is systematically lower for most 
sizes and in better agreement with existing theory [201 128] . 
N — 8 and N = 18 are minima as expected for a closed 
shell. There is also a deep minimum at TV = 10, (although 
the variation between experimental runs is larger for TV = 10). 
This supports the prolate structure for TV = 10 as predicted 
from the Clemenger-Nilsson model (lower) a/N for Najv 
(N = 10 — 200) The shell closings have been marked. They 
coincide with the extrema of the oscillations about the de- 
scending trend. At N = 200 the clusters are still far from the 
polarizability of bulk Na metal which is 9.4 A 3 /TV. Shown for 
comparison is the prediction of the spherical jellium model 
due to Ekardt g] 

These oscillations in the polarizability with the shell 
structure were predicted by Ekardt for Na |4j, and by 
Puska and co-workers for Li and Al clusters [51 [5] in 
the spherical jellium approximation. Puska et. al. [5] 
qualitatively explain this behavior as follows. By defini- 
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tion, for an electron in a quantum state with principal 
quantum number 1, there are no other electrons with 
the same angular momentum and lower principal quan- 
tum number. Consequently, electrons in these shells do 
not experience the Pauli repulsion from electrons in pre- 
viously occupied shells with the same angular momen- 
tum. Therefore, the orbitals of these electrons penetrate 
deeper into the cluster and their spillout is reduced. In 
contrast, electrons in shells with principal quantum num- 
ber 2 experience the Pauli repulsion from electrons with 
identical angular momentum in a previously filled shell 
(for example, electrons in the 2d shell are repelled by 
electrons in the Id shell.) This repulsion enhances the 
spillout and causes a/N to increase as this shell is filled. 

A triaxial distortion can also enhance the axis aver- 
aged a of a cluster. However, estimates of the magni- 
tude of this effect using values of the distortion parame- 
ter from photoabsorption experiments [7 shows that it is 
too small to account for the magnitude of the oscillations. 
It is also noteworthy that, a significant anomaly in the 
generally smooth trend is observed at TV = 55. This sud- 
den drop in a/N immediately before the electronic shell 
closing at N — 58, and is likely caused by the geometric 
shell closing. 

Overall the polarizabilities are systematically smaller 
than reported in previous experiments [11H13) . and are 
in closer agreement with existing theory 20, 28 . This 
effect has been predicted [2(5] and is related to thermal 
expansion. It is surprising that theory gives a to within 
5-10%, while the dipole moments are off by orders of 
magnitude, but this was already nearly the case with the 
spherical jellium model which has zero dipole moment by 
symmetry 0]. 

In conclusion, the electric deflection measurement dis- 
cussed here gives a comprehensive picture of the response 
of small sodium clusters to static electric fields. The 
nearly vanishing electric dipole moments, even for clus- 
ters as small as the sodium trimer, demonstrates that 
the electric fields surrounding alkali clusters are very 
small, as expected for a classical metallic object. The ob- 
served dipole moments are much smaller than predicted 
by quantum chemical methods, indicating a fundamental 
challenge for the theoretical treatment of dipole moments 
in metallic clusters. 

We gratefully acknowledge helpful discussions with Xi- 
aoshan Xu. Funding was provided by the NSF under fund 
No. R7418-DMR-0605894 



[2] W. A. de Heer, Rev. Mod. Phys. 65, 611 (1993). 
[3] N. Lang and W. Kohn, Phys. Rev. B 1, 4555 (1970). 
[4] W. Ekardt et al., Solid State Communications 57, 661 
(1986). 

[5] M. Manninen, R. Nieminen, and M. Puska, Phys. Rev. 

B 33, 4289 (1986). 
[6] M. J. Puska, R. M. Nieminen, and M. Manninen, Phys. 

Rev. B 31, 3486 (1985). 
[7] M. Schmidt and H. Haberland, Eur. Phys. J. D 6, 109 

(1999). 

[8] H. Haberland et al, Phys. Rev. Lett. 94, 35701 (2005). 

[9] C. Hock et al, Phys. Rev. Lett. 102, 43401 (2009). 
[10] O. Kostko et al, Phys. Rev. Lett. 98, 43401 (2007). 
[11] W. D. Knight et al, Phys. Rev. B 31, 2539 (1985). 
[12] G. Tikhonov et al, Phys. Rev. A 64, 63202 (2001). 
[13] D. Rayane et al., Eur. Phys. J. D 9, 243 (1999). 
[14] R. Moro et al., Science 300, 1265 (2003). 
[15] S. Yin et al., J. Supercond. Nov. Mag. 21, 265 (2008). 
[16] S. Schafer et al, J. Phys. Chem. A 112, 12312 (2008). 
[17] S. Schafer et al, J. Chem. Phys. 129, 044304 (2008). 
[18] V. Senz et al., Phys. Rev. Lett. 102, 138303 (2009), ISSN 
1079-7114. 

[19] Dipole moments were observed in Nb clusters (and 10 
alloys) up to N = 100. A rigorous theoretical explanation 
that can account for the loss of screening and all of the 
related experimental phenomena is currently lacking. See 
Ref. [TS] for more details and references to theoretical 
work. 

[20] I. A. Solovyov, A. V. Solovyov, and W. Greiner, Phys. 

Rev. A 65, 53203 (2002). 
[21] M. Broyer et al., C. R. Phys. 3, 301 (2002). 
[22] G. Bertsch, N. Onishi, and K. Yabana, Z. Phys. D Atoms 

34, 213 (1995). 

[23] M. A. El Rahim et al., J. Phys. Chem. A 109, 8507 
(2005). 

[24] J. Bowlan(2010), to be published. 

[25] We have applied corrections for other causes of broad- 
ening unrelated to a dipole moment: The dispersion in 
the beam velocity is negligible. The inhomogeneity of the 
deflection field causes a small but detectible broadening. 
This correction accounts for all of the broadening ob- 
served in the atomic beams used for calibration (e.g. Li, 
Al, In). 

[26] L. Coudert, W. Ernst, and O. Golonzka, J. Chem. Phys. 

117, 7102 (2002). 
[27] J. Paier, M. Marsman, and G. Kresse, J Chem. Phys. 

127, 024103 (2007). 
[28] S. Kiimmel et al., Eur. Phys. J. D 11, 239 (2000). 
[29] S. A. Blundell et al., Phys. Rev. Lett. 84, 4826 (2000). 



[1] J. Jackson, Classical Electrodynamics 3rd ed. (Wiley, 
1998). 



